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Abstract. Relativistic jets, formed in the vicinity of central supermassive black holes in AGN, show ample evidence 
connecting them to physical conditions in the accretion disc and broad-line region. The jets are responsible for a large 
fraction of non-thermal continuum emission (particularly during powerful flares), which makes understanding their 
physics an important aspect of studies of blazars characterised by profound flaring activity arising from extremely 
compact regions. Imaging and polarimetry of radio emission on milliarcsecond scales provided by very long baseline 
interferometry (VLBI) offers a range of possibilities for studying ultra-compact regions in relativistic jets and 
relating them to main manifestations of the blazar activity in AGN. Simultaneous monitoring of optical/high 
energy variability and evolution of parsec-scale radio structures yields arguably the most detailed picture of the 
relation between acceleration and propagation of relativistic flows and non-thermal continuum generation in blazars. 
Opacity effects provide a measure of magnetic field strength on scales down to ~ 1000 gravitational radii and trace 
the distribution of broad-line emitting material. Correlations observed between parsec-scale radio emission and 
optical and gamma-ray continuum indicate that a significant fraction of non-thermal continuum may be produced 
(particularly during flares) in extended regions of relativistic jet at distances up to 10 parsecs from the central 
engine. Combined with studies of jet component ejections and X-ray variability, these correlations also suggest 
that time delays, nuclear opacity, and jet acceleration may have a pronounced effect on the observed broad-band 
variability and instantaneous spectral energy distribution (SED). These effects will be reviewed below and discussed 
in the context of deriving accurate and self-consistent models for central regions of blazars. 



1. Introduction 

In the current astrophysical paradigm for active galactic 
nuclei (AGN), each constituent of an AGN contributes to 
a specific domain in the broad-band spectral energy distri- 
bution (Ghisellini & Tavecchio I2009|) . Variability of con- 
tinuum flux in AGN, signalling the activity of the central 
engine, is detected throughout the entire electromagnetic 
spectrum, on time-scales from days to years. Substantial 
progress achieved during the past decade in studies of ac- 
tive galactic nuclei (see Lobanov & Zensus 120061 for a re- 
view of recent results) has brought an increasingly wider 
recognition of the ubiquity of relativistic outflows (jets) in 
AGN (Falcke l200Tl Zensus [1997)1 and in blazars in partic- 
ular. 

Understanding the physics of blazars jets has acquired 
a particular importance after the launch of Fermi Gamma- 
Ray Space Telescope, as the compact, relativistic flows 
contribute strongly to the broad-band continuum - a fact 
that is still somewhat overlooked in physical models de- 
scribing high-energy emission from AGN. 

Emission properties, dynamics, and evolution of an 
extragalactic jet are intimately connected to the char- 
acteristics of the supermassive black hole, accretion disk 
and broad-line region in the nucleus of the host galaxy 
(Lobanov 2008) . The jet continuum emission is dominated 
by non-thermal synchrotron and inverse- Compton radi- 
ation (Unwin et al. 1997). The synchrotron mechanism 
plays a more prominent role in the radio domain, and 
the properties of the emitting material can be assessed 
using the turnover point in the synchrotron spectrum 



(Lobanov I1998b[ Lobanov & Zensus 1999jl , synchrotron 
self- absorption (Lobanov 1998a), and free-free absorption 
in the ambient plasma (Walker et al. 2000, Kadler et al. 

mm . 

There is growing evidence for relativistic flows con- 
tributing substantially to generation of non-thermal con- 
tinuum in the optical (Arshakian et al. 120101 Leon-Tavares 
et al. [20101 Jorstad et al. l2TJTU|) . X-ray ( Unwin et al. lTWl 
D'Arcangelo et al. [20071 Marscher et al. |2008)| Soldi et al. 
[MM Chatterjee et al.HEH), 7-ray (Otterbein et al. [15551 
Jorstad et al. 120101 Marscher et al. 120101 Savolainen et 
al. I2TJTU1 Schinzel et al. |2TJTU|) and TeV (Piner & Edwards 
I2001 Chariot et al. IMl Acciari et al. IMK?)) domains. 
Accurate spatial localisation of production sites of vari- 
able non-thermal continuum emission in AGN is therefore 
instrumental for understanding the mechanism for release 
and transport of energy in active galaxies. 

In the radio regime, very long baseline interferome- 
try (VLBI) enables direct imaging of spatial scales com- 
parable the gravitational radius, R g = G M^/c 2 , of the 
central black hole in AGN using ground VLBI observa- 
tions at 86 GHz and higher (c/., GMVA0 observations; 
Krichbaum et al. 2008) and space VLBI observations at 
centimetre wavelengths (Takahashi et al. l2004[) . Such high- 
resolution radio observations also access directly the re- 
gions where the jets are formed ( Junor et al. 1199 9;) , and 
trace their evolution and interaction with the nuclear en- 
vironment (Lobanov HZOOTl [2TJU51 and Middelberg & Bach 



Global 



Millimeter 



VLBI 



Array; 



http: / /www.mpifr-bonn. mpg.de/div / vlbi / globalmm 



1 



2 



Andrei Lobanov: Physical properties of blazar jets from VLBI observations 



2008). Evolution of compact radio emission from several 
hundreds of extragalactic relativistic jets is now system- 
atically studied with dedicated monitoring programs and 
large surveys using very long baseline interferometry (such 
as the 15 GHz VLBA@ survey (Kellermann et al. I2H04) 
the MOJAVeH survey (Lister & Homan [2005)) and a dedi- 
cated 22/43/86 GHz VLB A gamma-ray blazar monitoring 
program at the Boston UniversitjQ (Jorstad et al. I200ip . 
These studies, combined with optical and X-ray studies, 
yield arguably the most detailed picture of the vicinity of 
supermassive black holes in AGN (Marscher 2005). 

Presented below is a brief (and certainly incomplete) 
summary of recent results from VLBI studies of compact 
extragalactic radio sources, outlining the physical proper- 
ties of relativistic parsec-scale jets and their relation to 
supermassive black holes, accretion disks and broad-line 
regions in prominent blazars. 

2. Physics of compact jets 

Jets in active galaxies are formed in the immediate vicinity 
of the central black hole (Camenzind 2005), at distances of 
10-10 2 R g (Meier |2"009"|) . The jets carry away a fraction of 
the angular momentum and energy stored in the accretion 
flow (Blandford & Payne IT9821 Hujeirat et al. [20031 or 
corona (in low luminosity AGN; Merloni & Fabian [2002p 
and in the rotation of the central black hole (Blandford & 
Znajek 1977, Koide et al. 120021 Komissarov l2"001)l Semenov 
et al. 12002). 

The production of highly-relativistic outflows requires 
a large fraction of the energy to be converted to Poynting 
flux in the very central region (Sikora et al. I2005p . The 
efficiency of this process may depend on the spin of the 
central black hole (Meier I1999[) . The collimation of such 
a jet requires either a large scale poloidal magnetic field 
threading the disk (Spruit et al. U997[) or a slower and more 
massive MHD outflow launched at larger disk radii by cen- 
trifugal forces (Bogovalov & Tsinganos 2005). The flowing 
plasma is likely to be dominated by electron-positron pairs 
(Wardle et al. [TOOSl Hirotani |2"0P5|) although a dynami- 
cally significant proton component cannot be completely 
ruled out at the moment (Celotti & Fabian [1993 ). 

Acceleration or collimation of the flow may be com- 
plete within about 10 3 R g (Meier et al. I2009P or continue 
all the way to scales of a few parsecs (Vlahakis & Konigl 
2004]). At distances of ~ 10 3 i? g , the jets become visible 
in the radio regime. Recent studies indicate that at 10 3 - 
10 5 R g (< 1 pc) the jets are likely to be dominated by pure 
electromagnetic processes such as Poynting flux (Sikora 
et al. I2005P or have both MHD (kinetic flux) and electro- 
dynamic components (Meier [5003) . At larger scales, the 
jets are believed to be kinetic flux-dominated. The mag- 
netic field is believed to play an important role in accel- 
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erating and collimating extragalactic jets on parsec scales 
(Vlahakis & Konigl 2004}. Three distinct regions with dif- 
ferent physical mechanisms dominating the observed prop- 
erties of the jets can be considered: 1) ultra-compact jets 
(on scales of up to ~ 1 pc) where collimation and accel- 
eration of the flow occurs; 2) parsec-scale flows (^10pc 
scales) dominated by relativistic shocks and 3) larger-scale 
jets (~ 100 pc) where plasma instability gradually becomes 
dominant. This picture may be further complicated by 
transverse stratification of the flow, with the jet velocity, 
particle density and magnetic field changing substantially 
from the jet axis to its outer layers. As a practical result 
of this stratification, shocks and instabilities may in fact 
co-exist on all scales in the jets, with instabilities simply 
remaining undetected in compact flows, owing to limited 
resolution and dynamic range of VLBI observations. 

2.1. Ultra-compact jets 

Blazar jets usually feature a bright, compact (often un- 
resolved) "core" and a weaker, extended jet (often also 
transversely unresolved), with several regions of enhanced 
emission (traditionally branded "jet components" ) embed- 
ded in the flow and separating from the core at appar- 
ently superluminal speeds. The radio core is offset from 
the true base of the jet, and this "invisible" flow is proba- 
bly Poynting flux-dominated (Meier [TOOOl |2"0031 Sikora et 
al. 121705)1 . 

The radio core has a flat spectrum, expected to re- 
sult from synchrotron self-absorption in a conically ex- 
panding ultra-compact flow (Konigl 1981). As a result, 
the observed position, r c , of the core depends on the fre- 
quency of observation, v, so that r c oc j/~ 1 / fcr (this is so 
called "core shift" effect). If the core is self-absorbed and 
in equipartition, the power index fc r = 1 (Blandford & 
Konigl 1979). Recent measurements (Kovalev et al. 2008) 
have also shown that the frequency dependent core shift 
increases during flares as expected from the synchrotron 
self-absorption. 

The core shift can be used for determining basic physi- 
cal properties of the ultra-compact flow and the surround- 
ing absorbing material (Lobanov |1998al) . Changes of the 
core position measured between three or more frequencies 
can be used for determining the value of k r and estimating 
the strength of the magnetic field, B corc , in the nuclear re- 
gion and the offset, i? C orc, of the observed core positions 
from the true base of the jet. The combination of B core 
and i?core gives an estimate for the mass of the central 
black hole M bh a 7 x 10 9 M o (B corc /G) 1 / 2 (i? corc /pc) 3 / 2 . 

Core shift measurements provide estimates of the total 
(kinetic + magnetic field) power, the synchrotron luminos- 
ity, and the maximum brightness temperature, lb, max in 
the jets can be made. The ratio of particle energy to mag- 
netic field energy can also be estimated, from the derived 
lb, max- This enables testing the original Konigl model and 
several of its later modifications (e.g., Hutter & Mufson 
1986; Bloom & Marscher 1996). The known distance from 
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the nucleus to the jet origin can also enable constrain- 
ing the self-similar jet model (Marscher 1995) and the 
particle-cascade model (Blandford & Levinson [T995 ) . 

Studies of free-free absorption in AGN indicate the 
presence of dense, ionised circumnuclear material with 
T c « 10 4 K distributed within a fraction of a parsec of 
the central nucleus (Lobanov I1998al Walker et al. I2000p . 
Properties of the circumnuclear material can also be stud- 
ied using the variability of the power index k r with fre- 
quency. This variability results from pressure and den- 
sity gradients or absorption in the surrounding medium 
most likely associated with the broad-line region (BLR). 
Changes of k r with frequency, if measured with required 
precision, can be used to estimate the size, particle den- 
sity and temperature of the absorbing material surround- 
ing the jets. Estimates of the black hole mass and size of 
the BLR obtained from the core shift measurements can 
be compared with the respective estimates obtained from 
the reverberation mapping and applications of the Mbh - 
a* relation. 

The non-thermal continuum radio emission from the 
jet core does not appear to have strong shocks (Lobanov 
1998b), and its evolution and variability can be explained 
by smooth changes in particle density of the flowing 
plasma, associated with the nuclear flares in the central 
engine (Lobanov & Zensus IT999p . Compelling evidence ex- 
ists for acceleration (Vlahakis & Konigl [2004[ Bach et al. 
[20051 Lee et al. [2008) and collimation ( Junor et al. [T9991 
Krichbaum et al. 2008) in the ultra-compact flows. 

The brightness temperature of the radio emission from 
the cores reaches the inverse- Compton limit of w 10 12 K 
(Unwin et al. 119971 Lobanov et al. 120001 Kovalev et al. 
2005), while it drops rapidly to the equipartition limit 
of rs 5 x 10 10 K in the jet components moving down- 
stream from the core (Lahteenmaki et al. [1999, Lobanov 
et al. 120001 Homan et al. 12006|) . This supports earlier con- 
clusions that ultra-compact jets are particle-dominated, 
while the plasma in moving jet components is likely to 
be close to the equipartition (Unwin et al. 119971 Lobanov 
1998aJ Hirotani). Combining these calculations with esti- 
mates of the jet kinetic power provides strong indications 
that the relativistic fraction of the outflowing material is 
most likely represented by the electron-positron plasma 
(Reynolds et al. 115951 Hirotani et al. [20001 Hirotani [5003) • 

2.2. Parsec-scale flows: shocks and instabilities 

Parsec-scale outflows are characterised by pronounced 
curvature of trajectories of superluminal components 
(Kellermann 120041 Lobanov & Zensus I1999[) . rapid 
changes of velocity and flux density (Lister et al. l2009[) and 
predominantly transverse magnetic field 120051 Statistical 
studies of speed and brightness temperature distributions 
observed in the superluminal features propagating on par- 
sec scales indicate that the jet population has an envelope 
Lorentz factor of « 30 and an un-beamed luminosity of 
- 1 x 10 25 WHz^ 1 (Cohen et al. 1200711 . 



Physical conditions of the jet plasma can be assessed 
effectively by studying the spectral peak (turnover point) 
of the synchrotron emission. Mapping the turnover fre- 
quency distribution provides a sensitive diagnostic of the 
plasma (Lobanov I1998b|) . Obviously, the observed mor- 
phology and velocity of these flows are affected substan- 
tially by Doppler boosting, aberration, and time delays, 
which makes uncovering true physical properties a non- 
trivial task (Gomez et al. [T0M1 Cohen et al. [20071) . As 
a result, it becomes difficult to distinguish between ap- 
parent and true physical accelerations of the moving fea- 
tures (Lister et al. I2009[) and making such a distinction 
often requires a detailed physical modelling of a given jet 
component (Lobanov & Zensus 1999, Homan et al. I2003p . 
Similarly to stellar jets, rotation of the flow is expected 
to be important for extragalactic jets (Fendt I1997p . but 
observational evidence remains very limited on this issue. 

VLBI studies have demonstrated that relativistic 
shocks are prominent in jets on parsec scales, which is 
manifested by strong polarization (Ros et al. 2000) and 
rapid evolution of the turnover frequency of synchrotron 
emission (Lobanov et al. 11997} Lobanov & Zensus 1999). 
Evidence is growing for the presence of stationary features 
in parsec-scale flows, typically separated by ~ 1 pc dis- 
tance from the jet core (Kellermann et al. 120041 Savolainen 
et al. [20001 Lister et al. [20051 Arshakian et al. l20T01 Leon- 
Tavares et al. l2PTP|) . 

Specific geometric conditions and extremely small 
viewing angles could lead to formation of stationary fea- 
tures in relativistic flows (Alberdi et al. I2000|) . However, 
a more general and physically plausible explanation is 
offered by standing shocks (for instance, recollimation 
shocks in an initially over-pressurised outflow; Daly & 
Marscher [TQggl Gomez et al. [TOOSI Perucho & Marti 
2007). Such standing shocks may play a major role in ac- 
celerating particles near the base of the jet (Mandal & 
Chakrabarti 120081 Becker et al. I2008p . and could be re- 
sponsible for the persistent high levels of polarization in 
blazars (D'Arcangelo et al. [20071 Marscher et al. |2008)| ). 
More speculatively, the stationary features in jets could 
also be the sites of continuum emission release due to 
conversion from Poynting flux-dominated to kinetic flux- 
dominated flow. 

Complex evolution of the moving shocked regions is 
revealed in observations (Gomez et al. 120011 Jorstad et al. 
2005, Lobanov & Zensus IT999P and numerical simulations 
(Agudo et al. l200ip of parsec-scale outflows. However, the 
shocks are shown to dissipate rapidly (Lobanov & Zensus 
1999), and shock dominated regions are not likely to ex- 
tend beyond ~ 100 pc. Starting from these scales, instabil- 
ities (most importantly, Kelvin- Helmholtz instability, cf., 
Hardee 2000) determine at large the observed structure 
and dynamics of extragalactic jets (Lobanov et al. [1998, 
Walker et al. [20011 Lobanov & Zensus [20011 Lobanov et 
al. [20031 Hardee et al. [20051 Perucho et al. 12005)1 . 

The elliptical mode of the instability is responsible for 
appearance of thread-like features in the jet interior , while 
overall oscillations of the jet ridge line are caused by the 
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helical surface mode (Lobanov & Zensus 2001). Successful 
attempts have been made to represent the observed bright- 
ness distribution of radio emission on these scales, using 
linear perturbation theory of Kelvin-Helmholtz instabil- 
ity (Lobanov & Zensus 120011 Lobanov et al. 2003] Hardee 
et al. I2005|) and a spine-sheath (analogous to the two- 
fluid) description of relativistic flows (Canvin et al. 2005, 
Laing & Bridle 2004). Non-linear evolution of the insta- 
bility (Perucho et al. I2004al I2004b|) . stratification of the 
flow (Perucho et al. 120051) . and stabilisation of the flow 
via magnetic field (Hardee 120071) are important for repro- 
ducing the observed properties of jets. At larger scales, 
the helical surface mode of Kelvin-Helmholtz instability is 
likely to be one of the most important factor for disrupt- 
ing and destroying the outflows (Lobanov et aL I2U031 l^006, 
Perucho & Lobanov [2"007[) . 

2.3. Structure of the magnetic field 

The structure of magnetic field in blazar jets can be as- 
sessed with VLBI via linear (e.g., Lister & Homan 2005) 
and circular polarisation (Homan & Lister [2006) measure- 
ments. The ultra-compact jets (VLBI cores) are shown to 
be typically less than 5 % linearly polarised, with the po- 
larisation angle in BL Lac objects showing a stronger ten- 
dency to be aligned with the inner jet direction (Lister & 
Homan 2005). The low degree of polarisation in the VLBI 
cores can be caused by a disordered magnetic field (Hughes 
2005} or strong Faraday de-polarisation (Zavala & Taylor 
2004), or result from the "beam de-polarisation" if the 
magnetic field in the cores is structured on scales much 
smaller than the resolution of the VLBI experiments. 

The fractional linear polarisation in moving jet compo- 
nents grows with increasing distance from the core (Lister 
& Homan I2005|) , and the position angle of the polarisa- 
tion vector is again better aligned with the jet direction 
in the BL Lac objects. At the same time, the linearly 
polarised emission in transversely resolved flows displays 
both strong rotation of the polarisation angle near the core 
and remarkable edge brightening (e.g., Ros et al. 2000), 
with the polarisation vectors tending to be perpendicular 
to the jet direction. This indicates overall complexity and 
likely stratification of the magnetic field in the jets, with 
the internal part ("spine") of the flow being dominated 
by a helical magnetic field further compressed by rela- 
tivistic shocks, while the external, slower moving layers 
( "sheath" ) of the flow could be dominated by a longitudi- 
nal magnetic field. This conclusion is in a good agreement 
with the turnover frequency distribution (Lobanov et al. 
I1997p . internal structure (Lobanov & Zensus [200ip . and 
transverse velocity stratification in the jets (Perucho et al. 
2006). Here again, the effect that the flow, and magnetic 
field, rotation (due to residual angular momentum inher- 
ited from the initial disk-to-jet coupling) remains difficult 
to assess. Presence of a strong toroidal or helical magnetic 
field either in the spine or in the sheath has also been 
suggested from observations of Faraday rotation gradients 



across the flow (Asada et al. [20021 Zavala & Tavlor [20051 
Attridge et al. [277051 Gomez et al.[Mp. 

Circular polarisation has been detected in a number 
of AGN jets (Homan & Lister I2006[) , with a typical level 
of polarisation of <0.5%. The circular polarisation can 
be either intrinsic to the synchrotron emission (Legg & 
Westfold 119681 implying the presence of a strong rela- 
tivistic proton component) or result from scintillations 
(Macquart & Melrose 2000), relativistic effects in disper- 
sive plasma (Broderick & Blandford 2002) or Faraday con- 
version of linearly polarised synchrotron emission from 
electron-positron plasma (Jones & Q'Dell lT977p . The ob- 
served properties of circularly polarised emission in blazar 
jets support the last mechanism for its formation (Wardle 
et al.[nSBl Homan & Lister 120061 Homan et al. 1MJ9")) . 

The strength of the magnetic field in blazar jets 
is typically assessed by combining multi-band measure- 
ments (Unwin et al. 1997), by deriving information about 
the peak in the synchrotron spectrum (Marscher 1983, 
Lobanov et al. 119971 Lobanov I1998bi Savolainen et al. 
2008, Sokolovsky et al. I2010|) . or by using the opacity due 
to synchrotron self-absorption (Lobanov 1998a) . In most 
objects, values around 1 G are obtained for the VLBI cores 
and lower magnetic field is measured in the jets, all falling 
well in agreement with magnetic field generation being ul- 
timately processes in the magnetised accretion disk (Field 
& Rogers H993[ ). 

2.4. Periodic changes of the structure 

Structural changes are abound on milliarcsecond scales in 
the blazar jets, enhanced and magnified by small viewing 
angles and relativistic effects. In addition to extreme cur- 
vature observed in the jet ridge line of some objects (e.g., 
Polatidis et al. U995i , the position angle of the inner jet 
(as traced by jet components nearest to the core) changes 
substantially both as a function of observing frequency 
(Savolainen et al. 120061 Agudo et al. I2007p and in time 
(Mutel & DennCSJOSl Lobanov & Roland [HJS) . 

The frequency dependent changes are most likely 
caused by the opacity and spectral index gradients in the 
flow. The temporal variations of the position angle can re- 
sult from precession and rotation of the flow (Camenzind 
et al. I1992p as well as from the pattern motion of Kelvin- 
Helmholtz instability (Hardee [2UTJ51 Hardee et al. |2"0"0"5[) . 
in which case the motion should also be evident in long- 
term evolution of the ridge line of a flow (e.g., Krichbaum 
et al. l200T[) . 

Variations of the jet position angle, as well as the ob- 
served morphology of parsec-scale jets and trajectories of 
superluminal features propagating in the jets, are often 
described in terms of a helical geometry (Steffen et al. 
H), with helicity supposed to be arising from some peri- 
odic process in the nucleus. Jet precession, both in sin- 
gle and binary black hole systems, have been commonly 
sought to be responsible for the observed helicity on par- 
sec scales. However, as the observed periods of the posi- 
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tion angle changes (and, similarly, periods inferred from 
the component trajectories and oscillations of the jet ridge 
lines) are typically within a range of a few years, the 
precession-based models face severe difficulties (Lobanov 
& Roland 2005) as they require either allowing for ex- 
tremely small (< 10 4 i? g ) orbital separations in supermas- 
sive binary black holes or adopting an assumption that 
the jet direction responds exclusively to changes in the 
innermost parts of the accretion disk and it is decoupled 
from the spin of the central black hole. In view of these 
difficulties, rotation of the flow and pattern motion of the 
instability seem to be more viable alternatives. The preces- 
sion of the flows, evident on kiloparsec-scales (Gower et al. 
1982, Hardee et al. I1994p acts on much longer timescales 
(typically > 10 4 yr), typically of several hundred years 
and longer, and thus should only be visible in long-term 
changes of the position angle of the entire milliarcsecond- 
scale jet (Lobanov & Roland I2005p . 

3. Emission from blazar jets 

VLBI observations enable tracing both temporal and spa- 
tial changes in radio emission from blazar jets, offering a 
unique opportunity for connecting these changes to prop- 
erties of the blazar emission observed in other wavelength 
domains and even localising spatially the dominant com- 
ponents in the broad-band emission. Emission from jets 
may also contribute substantially to the broad-band SED 
(Yuan et al. I2002p . alongside with the canonical contri- 
butions from the accretion disk and hot material in the 
vicinity of the central black hole. VLBI observations help 
here enormously by providing an accurate measure of ra- 
dio emission produced in these regions and excluding con- 
tributions from kiloparsec-scale jets and lobes of radio 
sources. 

3.1. Parsec-scale radio emission 

Parsec-scale radio emission is variable on timescales from 
decades (e.g., O'Dea et al. 119841 Asada et al. I2U06|) to 
hours (Savolainen & Kovalev I2008p . with the longest 
timescales most likely related to large-scale changes in 
the nuclear region feeding the black hole, and the short- 
est timescales resulting from interstellar scintillations. The 
variability on hour-to-day timescales may also be related 
to "quasar QPO" type of variations recently observed in 
the soft X-ray band (Gierlinski et al. I2008|) . remembering 
that if variability periods scale with the black hole mass, 
this would correspond to canonical QPO seen in X-ray 
binaries at frequencies around 100 Hz. 

On intermediate timescales (months-to-years), most 
of variable radio emission is believed to be associated 
with flares in the VLBI cores (Lobanov & Zensus 1999) 
and shock evolution of plasma propagating downstream 
(Hughes et al. [1985, Marscher & Gear 119851 Marscher 
1990). The radio flares last, on average, for 2.5 years 
(at a wavelength of about 1cm; Hovatta et al. I2008p . 
and in many objects they are repeated quasi-periodically 



(Hovatta et al. 2007). The flares are firmly associated with 
ejections of new jet components (Valtaoja et al. 119991 
Lobanov & Zensus [M)9l Marscher et al. 120021 Chatterjee 
et al. I2009p . but the release of non-thermal continuum 
emission may not necessarily be restricted to the vicin- 
ity of the black hole or even the radio core (Arshakian et 
al. [20101 Leon-Tavares et al. I20T01 Schinzel et al. [2010)1 . 
Detailed evolution of a flaring emission is best determined 
through variations of its turnover peak (Otterbein et al. 

1998, Fromm et al. I2010p and can be further constrained 
using observed kinematic properties of an emitting region 
(Lobanov & Zensus [T999|) . 

Quasi-periodic variability of the radio emission from 
the ultra-compact jets is most likely related to instabilities 
and non-stationary processes in the accretion disks around 
central black holes in AGN (Igumenschev & Abramowicz 

1999, Lobanov & Roland 2005, Alternative explanations 
involve binary black hole systems in which flares are 
caused by passages of the secondary through the accre- 
tion disk around the primary Ivanov et al. 1998, Lehto & 
Valtonen 119961 Similarly to the attempts of using binary 
black holes to explain short-term morphological changes, 
these models require very tight binary systems, with or- 
bits of the secondary lying well within 10 4 R g of the pri- 
mary, which poses inevitable problems for maintaining 
an accretion disk around the primary (for massive secon- 
daries; Lobanov 2008) or rapid alignment of the secondary 
with the plane of the accretion disk (for small secondaries; 
Ivanov et al. UM^ . 

3.2. Blazar jets and broad-band continuum 

Relativistic flows are prominent emitters in all bands of 
the electromagnetic spectrum, generating optical and X- 
ray emission even on kiloparsec scales (Schwartz et al. 
[217001 Marshall et al. [2UD21 Siemiginovska et al. |2"U0"21 
Sambruna et al. I2008P and at TeV energies (Acciari et al. 
I2010p . The jet plasma is believed to emit via synchrotron 
emission in the radio to soft X-ray range and via inverse 
Compton emission in the hard X-ray to TeV range (Acciari 
et al. 120101 Marscher et al. I2010p . Contested is, however, 
the primary source of the seed photons for the inverse 
Compton part of the radiation. These could be the syn- 
chrotron photons themselves (synchrotron self Compton 
mechanism, SSC) or photons from an external radiation 
field located, for instance, in the accretion disc emission, 
X-ray corona, the broad-line region, the infrared emitting 
torus and the cosmic background radiation (Ghisellini & 
Tavecchio [20091) . 

Results from the first year of the Fermi /L AT opera- 
tions have enabled most detailed studies of the connection 
between relativistic flows and 7-ray production in blazars. 
Early statistical comparisons of the properties of 7-ray 
emission and compact radio jets in blazars indicate un- 
equivocally that they are closely related (Pushkarev et 
al. 120091 Savolainen et al. I2010p . Several of the 7-ray 
flares detected with Fermi /LAT can be associated with 
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emission from accelerated plasma cloud embedded in the 
jets (Marscher et al l20T01 Schinzel et al. |2"UTU|) . although 
detailed localisation of the 7-ray emitting sites remain 
elusive, clearly calling for continuation of extensive, co- 
ordinated Fermi /LAT and VLBI campaigns. 

There is now growing evidence for the broad-band con- 
tinuum (and its flaring components in particular) to be 
produced at multiple locations in AGN (Arshakian et al. 
2010), with the emission in different bands dominated 
by contributions from spatially different regions (Leon- 
Tavares 2010). These findings put an additional strain on 
single-zone models commonly used for fitting the broad- 
band SED in AGN. The situation with the quiescent (and 
slowly variable) component of the 7-ray emission is also 
puzzling, with indications that it may be produced in the 
region of the jet extending up to ~10pc (Schinzel et al. 

EHUD- 

These results give a compelling indication that forma- 
tion of jet components may correspond to the strongest 
nuclear events, while some of them may not survive a 
passage through a standing shock (Leon-Tavares et al. 
2010]). The source of the continuum emission is localised 
not only in the accretion disk (at the extreme vicinity of 
the black hole) but also in the entire acceleration zone of 
the jet, with strong flares happening both near the cen- 
tral black hole and at a standing shock in the jet. Clearly, 
a more general and systematic study of relation between 
the radio, optical and X-ray emission in blazars is strongly 
justified, and co-ordinated VLBI-Ferrm/LAT campaigns 
would one of the prime methods for such studies. 

4. Outlook 

More than thirty years after their appearance on the sci- 
entific scene, blazars remain one of the focal points of 
extragalactic astrophysics. High-resolution radio observa- 
tions of blazars provide essential information about struc- 
ture, kinematics and emission of relativistic flows in these 
objects on scales inaccessible to direct imaging in other 
bands. This information turns out to be arguably an in- 
dispensable tool for constructing viable physical models 
capable of explaining the blazar phenomenon and its spec- 
tacular observational manifestations. 

There are still a number of unanswered questions 
about the jets themselves and their relations to the broad- 
band emission produced in blazars. Extragalactic jets are 
an excellent laboratory for studying physics of relativistic 
outflows and probing conditions in the central regions of 
active galaxies. Recent studies of extragalactic jets show 
that they are formed in the immediate vicinity of central 
black holes in galaxies and carry away a substantial frac- 
tion of the angular momentum and energy stored in the 
accretion flow and rotation of the black hole. The jets are 
most likely collimated and accelerated by electromagnetic 
fields. Relativistic shocks are present in the flows on small 
scales, but dissipate on scales of decaparsecs. Plasma in- 
stabilities dominate the flows on larger scales. Convincing 
observational evidence exists, connecting ejections of ma- 



terial into the flow with instabilities in the inner accretion 
disks. 

In the coming years, new breakthroughs should be 
coming in studies of relativistic jets and blazars, enriched 
by Fermi /LAT results and by an effective use of the VLBI 
potential for studies of blazars. There are several potential 
focal points for these studies. 

In particular, the jet composition remains an open is- 
sue, in particular with regard to the role of relativistic 
protons in high-energy emission production in extreme 
vicinity of the central black holes (while it seems that 
pair plasma is chiefly responsible for the emission from 
parsec-scale jets). 

Accurate spatial localisation of the sites of high-energy 
continuum production will play a crucial role for modelling 
the broad-band SED of blazars and understanding their 
physical nature in general. To this end, variability of the 
high-energy continuum emission is best related to struc- 
tural and radiative changes in parsec-scale radio emission 
resolved by VLBI observations. 

Last but not least, the continued observational quest 
for reaching ever closer to the regions where the relativis- 
tic flows are formed will bring new answers not only about 
the physical nature of relativistic flows, but also about the 
physical properties of black holes and their connection to 
major manifestations of nuclear activity in galaxies. This 
would be particular important for studies of the 7-ray 
and TeV emission from blazars, as this emission has all 
chances to come from the immediate vicinity of the cen- 
tral black holes in galaxies. Combination of high-energy 
observations with VLBI extensive monitoring at centime- 
tre wavelengths and focused, specific VLBI programs at 
millimetre wavelengths certainly has a quality of the tool 
of choice for such studies. 
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